Elliptic Flow: Transition from out-of-plane to in-plane Emission in Au + Au Collisions 
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We have measured the proton elliptic flow excitation function for the Au + Au system spanning 
the beam energy range 2-8 AGeV. The excitation function shows a transition from negative to 
positive elliptic flow at a beam energy, Etr ~ 4 AGeV. Detailed comparisons with calculations from 
a relativistic Boltzmann-equation are presented. The comparisons suggest a softening of the nuclear 
equation of state (EOS) from a stiff form (K~ 380 MeV) at low beam energies (EBeam < 2 AGeV) 
to a softer form (K~ 210 MeV) at higher energies (EBeam > 4 AGeV ) where the calculated baryon 
density p ~ 4po- 
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For many years, the investigation of the nuclear equation of state (EOS) has stood out as one of the primary driving 
forces for heavy ion reaction studies (e.g. IQQ)- Measurements of collective motion and, in particular, the elliptic 
flow have been predicted to provide information crucial for establishing the parameters of the EOS ]3| |5j. Theoretical 
conjectures have also focused on the notion that a transition to the quark-gluon plasma (QGP) is associated with a 
"softest point" in the EOS where the pressure increase with temperature is much slower than the energy density [0. 
Such a softening of the EOS is predicted to start at quark-antiquark densities comparable to those in the ground-state 
of nuclear matter tt, and also at relatively low temperatures if the baryon density is driven significantly beyond its 



normal value po [|S|lO|. At energies of 1 < EBeam ^11 AGeV, collision-zone matter densities are expected up to 
p ~ 6 — 8po 0>[QJ- Such densities could very well result in conditions favorable to a softening of the EOS. Therefore, 
it is important to investigate currently available elliptic flow data [in this energy range] to search for new insights into 
the parameters of the EOS and for any indication of its softening. 

Elliptic flow reflects the anisotropy of transverse particle emission at midrapidity. For beam energies of 1-11 AGeV 
this anisotropy results from a strong competition between "squeeze-out" and "in-plane flow" |j|,[5],D ■ The magnitude 
and the sign of elliptic flow depend on two factors: (a) the pressure built up in the compression stage compared to 
the energy density, and (b) the passage time of the projectile and target spectators. The characteristic time for the 
development of expansion perpendicular to the reaction plane can be estimated as ~ R/c s , where the speed of sound 
c s = y/ dp/de, R is the nuclear radius, p is the pressure and e is the energy density. The passage time is ~ 2R/('-fo vq), 
where vo is the cm. spectator velocity. Thus the "squeeze-out" contribution should reflect the ratio Cs/70 i>o @] which 
is responsible for the essentially logarithmic dependence of elliptic flow on the beam energy for ~ 1 < Eb eam < H 
AGeV §. 

Recent calculations have made specific predictions for the beam energy dependence of elliptic flow for Au + Au 
collisions at 1-11 AGeV ||. They indicate a transition from negative to positive elliptic flow at a beam energy E tr , 
which has a marked sensitivity to the stiffness of the EOS. In addition, they suggest that a phase transition to the 
QGP should give a characteristic signature in the elliptic flow excitation function due to significant softening of the 
EOS. In this Letter we present an experimental elliptic flow excitation function for the Au + Au system to establish 
E tr and to search for any hints of a softening of the EOS. 

The measurements were performed at the Alternating Gradient Synchrotron (AGS) at the Brookhaven National 
Laboratory. Beams of 197 Au (EBeam — 2, 4, 6, and 8 AGeV) || were used to bombard a 197 Au target of thickness 
calculated for a 3% interaction probability. Typical beam intensities resulted in ~ 10 spills/min with ~ 10 3 particles 
per spill. Charged reaction products were detected with the E895 experimental setup which consists of a time 
projection chamber (TPC) [H and a multisampling ionization chamber (MUSIC) Jl|. The TPC which was located 
in the MPS magnet (typically at 1.0 Tesla) provided good acceptance and charge resolution for charged particles 
— 1 < Z < 6 at all four beam energies. However, unique mass resolution for Z — 1 particles was not achieved for 
all rigidities. The MUSIC device, positioned ~ 10 m downstream of the TPC, provided unique charge resolution for 
fragments with Z > 7 for the 2 and 4 AGeV beams. Data were taken with a trigger for minimum bias and also for 
a bias toward central and mid-central collisions. Results are presented here for protons measured in the TPC for 
mid-central collisions. 

We use the second Fourier coefficient V2 = (cos 20), to measure the elliptic flow or azimuthal asymmetry of the 
proton distributions at midrapidity (|y C m| < 0.1) ; 

dN 

— ~ [1 + 2vi cos(</>) + 2v 2 cos(20)] . (1) 
dq> 

Here, <j> represents the azimuthal angle of an emitted proton relative to the reaction plane. The Fourier coefficient 
(cos 20) = 0, > 0, and < for zero, positive, and negative elliptic flow respectively. Measurements of V\ will be 
presented and discussed in a forthcoming paper H . 
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Our analysis proceeds in two steps. First, we determine the reaction plane and its associated dispersion for each 
beam energy. Second, we generate azimuthal distributions with respect to this experimentally determined reaction 
plane and evaluate (cos2</>). The vector Qi = Y^j=a w (Vj) Pj/Pj ^ s usec ^ to determine the azimuthal angle, $ p i a „ e , 
of the reaction plane |jl5|| . Here, and yj represent, respectively, the transverse momentum and the rapidity of 

baryon j (Z< 2) in an event. The weight w(yj) is assigned the value <pt> ; where p x is the transverse momentum in 
the reaction plane. < p x > is obtained from the first pass of an iterative procedure. 

The dispersion of the reaction plane as well as biases associated with detector efficiencies plays a central role 
in flow analyses Jl6|,[l7],[l9|] . Consequently, in Fig. [I] we show representative distributions for the experimentally 
determined reaction-plane (Qpiane), and the associated relative reaction-plane distributions ($12). The distributions 
have been generated for a mid-central impact parameter, i.e. multiplicities between 0.5 and 0.75 M max . Here, M max 
is the multiplicity corresponding to the point in the charged particle multiplicity distribution where the height of 
the multiplicity distribution has fallen to half its plateau value |pcfl . It is estimated that this multiplicity range 
corresponds to an impact parameter range ^5 — 7 fm. The $12 distributions (cf. Fig. [I]) which are important 



for assessing the role of the reaction-plane dispersion, have been obtained via the subevent method 15 1. That is, 
reaction planes were determined for two subevents constructed from each event; $12 is the absolute value of the 
relative azimuthal angle between these two estimated reaction planes. The essentially flat reaction plane distributions 
shown in Fig. |l|a reflect rapidity and multiplicity-dependent azimuthal efficiency corrections, applied to take account 
of the detection inefficiencies of the TPC. These corrections were obtained by accumulating the laboratory azimuthal 
distribution of the particles (as a function of rapidity and multiplicity) for all events and then including the inverse 
of these distributions in the weights for the determination of the reaction plane. The distributions shown in Fig. |l|a 
confirm the absence of significant distortions which could influence the magnitude of the extracted elliptic flow. The 
relative reaction-plane distributions ($12) shown in Fig. [l]b indicate mean values which increase with the beam energy 
from < $12 > /2 ~ 17.0° at 2 AGeV to ~ 36.1° at 8 AGeV. This increase suggests a progressive deterioration in 
the resolution of the reaction plane with increasing beam energy; however a reasonable resolution is maintained over 
the entire energy range. The $12 distributions serve as the basis for correcting the extracted elliptic flow values as 
discussed below. 

In Fig. |^, we show observed (or cf) 1 ) azimuthal distributions, for protons. The distributions, shown for several 
rapidity bins, have been generated for the same mid-central impact parameter range (~ 5 — 7 fm) discussed above. 
Several characteristic features are exhibited in Fig. |^. For example, as one moves away from midrapidity, the 
distributions exhibit shapes commonly attributed to collective sidewards flow. That is, for y > 0, the distributions 
peak at 0°, and, for y < 0, they peak at ±180°. Fig. || also shows that these anisotropies decrease with increasing 
beam energy. 

The primary feature of the midrapidity distributions contrasts with those obtained at other rapidities. At 2 AGeV, 
two distinct peaks can be seen at —90° and +90°. These peaks indicate a clear signature for the "squeeze-out" of 



nuclear matter perpendicular to the reaction plane 116 20-22] or negative elliptic flow. By contrast, at 6 and 8 AGeV, 



the midrapidity distributions peak at 0°, and ±180°. This latter anisotropy pattern is expected for positive elliptic 
flow. Thus, Fig. |t provides clear evidence for negative elliptic flow at 2 AGeV, positive elliptic flow for 6 and 8 AGeV, 
and near zero flow for Eseam = 4 AGeV. 

In order to quantify the proton elliptic flow, it is necessary to suppress possible distortions arising from imperfect 
particle identification (PId) . It is relevant to reiterate here that unique separation of ir + and protons was not achieved 
for all rigidities. To suppress such ambiguity we applied the following procedure. First, we plot the observed Fourier 
coefficient (cos 2(f) ) vs. pt with p t thresholds which allow clean particle separation (p t ~ 1 GeV/c). We then extract 
the coefficients for the quadratic dependence of (cos 2(f)) on p t (see inset in Fig. [|). These quadratic fits are restricted 
by the requirement that (cos 20'} = for p t — 0. Second, we correct the proton p t distributions for possible 7r + 
contamination by way of a probabilistic PId. The latter probabilities were obtained by extrapolating the exponential 
tails of the proton and n + rigidity distributions into the regions of overlap. A weighted average (relative number of 
protons in a p t bin times the (cos2<//) for that bin) was then performed to obtain (cos 2(f)) for each beam energy. 
Subsequent to this evaluation, we then use the relative reaction plane distribution at each beam energy (cf. Fig. |l|) 
to obtain dispersion corrections for the extracted Fourier coefficients [ p^| , p^7pl| ] . 

The relationship between the (cos 2(f)) (obtained with the estimated reaction plane) and the Fourier coefficient 
(cos 2(f>) relative to the true reaction plane is: 

(cos 2(f)) = (cos 2(p) (cos 2A$) . (2) 

where (cos2A$) is the correction factor determined from the (cos $12) flffl. Following the prescription outlined in 
Rcf. jl7|, we find correction factors which range from 0.79 at 2 AGeV to 0.29 at 8 AGeV. The correction factors are 
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summarized along with ((cos $12)) in Table 1. 

The corrected elliptic flow values, (cos 2<fi), are represented by filled stars in Fig. ||. This excitation function clearly 
shows an evolution from negative to positive elliptic flow within the region 2 < Eseam $ 8 AGeV and points to an 
apparent transition energy E tr ~ 4 AGeV. The solid and dashed curves represent the results of model calculations 
described below. Since the value of E tr is predicted to be sensitive to the parameters of the EOS B, it is important 
to examine additional constraints on its value. The inset in Fig. || shows the corrected (cos 20) values as a function 
of p t for protons. The solid curves in the figure represent quadratic fits to the data (2 and 6 AGeV) which are in 
agreement with the predicted quadratic dependence of (cos 20) on p t QJT^] . Of greater significance is the fact that a 
comparison of the pt dependence of the elliptic flow for 2, 4, and 6 AGeV, provides further direct evidence that the 
sign of elliptic flow changes as the beam energy is increased from 2 to 6 AGeV. The essentially flat pt dependence 
shown for 4 AGeV is consistent with E tr ~ 4 AGeV. 

To interpret these data, extensive calculations have been made to constrain the parameters of the EOS in the context 
of a newly developed relativistic Boltzmann-equation model (BEM) ||^3) . The phenomenological relativistic Landau 
theory of quasiparticles |^4| serves as a basis for the model which has nucleon, pion, A and N* resonance degrees of 
freedom as well as momentum dependent forces. Calculations were performed for both a soft (K = 210 MeV), and a 
stiff (K = 380 MeV) EOS for the same rapidity and impact parameter selections applied to the data. 

The elliptic flow excitation functions (calculated for free protons) are compared to the experimental data in Fig. ||. 
The dashed and solid curves represent the results for a stiff and a soft EOS respectively. In addition to the data from 
the present experiment (filled stars), Fig. || also shows experimental results for Au + Au reactions at 1.15 A GeV j25]j 
(filled triangle) and 10.8 A GeV P6| (filled circle). The experimental data are compatible with the excitation function 
predicted for a stiff EOS at beam energies 1 < EBeam < 2 AGeV. By contrast, the data show good agreement with 
the predictions for a soft EOS for 4 < EBeam $11 AGeV. This pattern is consistent with a softening of the EOS in 
semicentral collisions of Au + Au at ~ 4 AGeV. The calculated densities at maximum compression for these energies 
are of the order of ~ 4 po for the stiff EOS. 

In summary, we have measured an elliptic flow excitation function for mid-central collisions of Au + Au at 2, 4, 6, 
and 8 AGeV. The excitation function exhibits a transition from negative to positive elliptic flow with E tr ~ 4 AGeV. 
Detailed comparisons of these elliptic flow data have been made with calculated results from a relativistic Boltzmann- 
equation calculation. Within the context of a simple parametrization of the EOS, the calculations suggest an evolution 
from a stiff EOS (K~ 380 MeV) at low beam energies (< 2 AGeV) to a softer EOS (K~ 210 MeV) at higher beam 
energies (4 < EBeam $11 AGeV). Such a softening of the EOS could result from a number of effects, the most 
intriguing of which is the possible onset of a nuclear phase change ||. On the other hand, it should be noted that 
transport models have failed to reproduce low energy "squeeze-out" data with a single incompressibility constant 
p2| . Thus, additional experimental signatures as well as calculations based on other models will be necessary to 
test the detailed implications of these results. Nevertheless, the results presented here, clearly show that elliptic flow 
measurements can provide an important constraint on the EOS of high density nuclear matter. 
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Beam Energy (AGeV) 
2 


< cosA$ 12 > 
0.753 


< cos2(A$) > 
0.79 


4 


0.563 


0.62 


6 


0.359 


0.42 


8 


0.244 


0.29 



Table 1: Dispersion correction factors for each beam energy. 
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FIG. 1. Experimentally determined (a) reaction-plane (Qpiane) distributions, and (b) the associated relative reaction-plane 
distributions for 2, 4, 6, and 8 AGeV Au + Au. The reaction plane distributions include efficiency corrections for the 

TPC (see text). 
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FIG. 2. Azimuthal distributions (with respect to the reconstructed reaction plane) for 2, 4, 6, and 8 AGeV Au + Au. 
Distributions are shown for (a) —0.7 < y cm < —0.5, (b) —0.5 < y cm < —0.3, (c) —0.1 < y cm < 0.1, (d) 0.3 < y cm < 0.5, and 
(e) 0.5 < y C m. < 0.7. The mid-rapidiy selections for 4 - 8 AGeV also include a transverse momentum selection as indicated. 
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FIG. 3. Elliptic flow excitation function for Au + Au. The filled symbols represent the experimental data as indicated. 
The dashed curve (open circles) and the solid curve (open squares) represent the calculated excitation functions for a soft and 
a stiff EOS (both with momentum dependence) respectively. The inset shows the [dispersion corrected] transverse momentum 
dependence of the elliptic flow for the 2, 4 and 6 AGeV beams. 
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